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Type la supernovae are thought to be thermonuclear explosions of accreting 
white dwarfs that reach a critical mass limit. Despite their importance as cos- 
mological distance indicators, the nature of their progenitors has remained 
controversial. Here we report the detection of circumstellar material in a nor- 
mal Type la supernova. The expansion velocities, densities and dimensions 
of the circumstellar envelope indicate that this material was ejected from the 
progenitor system. The relatively low expansion velocities appear to favor a 
progenitor system where a white dwarf accretes material from a companion 
star which is in the red-giant phase at the time of explosion. 

Due to their extreme luminosities and high homogeneity, Type la Supernovae have been 
used extensively as cosmological reference beacons to trace the evolution of the Universe ( 1, 
2). However, despite significant recent progress, the nature of the progenitor stars and the 
physics which govern these powerful explosions have remained very poorly understood (3, 4). 
In the presently favored single-degenerate model, the supernova progenitor is a white dwarf 
in a close binary accreting from a non-degenerate companion (5); the white dwarf explodes in 
a thermonuclear explosion when it approaches the Chandrasekhar limit. A direct method for 
investigating the nature of the progenitor systems of Type la supernovae (hereafter SNe la) is to 
search for signatures of the material transferred to the accreting white dwarf in the circumstellar 
material (CSM). Previous attempts have aimed at detecting the radiation which would arise 
from the interaction between the fast moving SN ejecta and the slow moving CSM in the form 
of narrow emission Unes (6), radio (7) and X-ray emission (8). The most stringent upper limit to 
the mass loss rate set by radio observations is as low as 3 x 10~^ solar masses per year (M© yr~^) 
for an assumed wind velocity of lOkms"^ (7). Two remarkable exceptions are represented by 
two peculiar SNe la, SN 2002ic and SN 2005 gj, which have shown extremely pronounced 
hydrogen emission lines (9, 10), that have been interpreted as a sign of strong ejecta-CSM 
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interaction (11). However, the classification of these supemovae as SNe la has recently been 
questioned (12), and even if they were SN la, they are unlikely to account for normal Type la 
explosions (7) that, so far, lack any signature of mass transfer from a hypothetical donor. Here, 
we report direct evidence of CSM in a SN la that has shown a normal behavior at X-ray, optical 
and radio wavelengths. 

SN 2006X was discovered in the Virgo Cluster spiral galaxy NGC 4321 (13). A few days 
later, the object was classified as a normal Type la event occurring 1-2 weeks before maximum 
light, which was affected by substantial extinction (14). Prompt Very Large Array (VLA) ob- 
servations have shown no radio source at the SN position (15), establishing one of the deepest 
and earliest Umits for radio emission from a Type la, and implying a mass-loss rate of less than 
a few 10~^ Mq yr~^ (for a low wind velocity of lOkms"^). The SN was not visible in the 
0.2-10 keV X-rays band down to the SWIFT sateUite detection limit (8). 

We have observed SN 2006X with the Ultraviolet and Visual Echelle Spectrograph mounted 
at the European Southern Observatory 8.2m Very Large Telescope. Observations were carried 
out on four different epochs, which correspond to days —2, +14, +61 and +121 with respect 
to S-band maximum Ught. Additionally, a fifth epoch (day +105) was covered with the High 
Resolution Echelle Spectrometer mounted at the 10m Keck telescope (16). The most remark- 
able finding from our data is the clear evolution seen in the profile of the Na I D lines (5889.95, 
5895.92 A). In fact, besides a strongly saturated and constant component, arising in the host 
galaxy disk (see section S2, Fig. SI), a number of features spanning a velocity range of about 
100 km s~^ appear to vary significantly with time (Fig. 1, Fig. S2). SN 2006X is projected onto 
the receding side of the galaxy, and the component of the rotation velocity along the line of sight 
at the apparent SN location is about +75 km s~^ (17), which coincides with the strongly satu- 
rated Na I D component, the saturated Ca II H&K lines, and a weakly saturated CN vibrational 
band (0-0) (Fig. 2, Fig. SI). This and the lack of time evolution proves that the deep absorption 
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Figure 1: Time evolution of the Sodium D2 component region as a function of elapsed time 
since 5-band maximum light. The heliocentric velocities have been corrected to the rest-frame 
using the host galaxy recession velocity. All spectra have been normalized to their continuum. 
In each panel, the dotted curve traces the atmospheric absorption spectrum. 
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arises within the disk of NGC 4321 in an interstellar molecular cloud (or system of clouds) that 
is responsible for the bulk of the reddening suffered by SN 2006X (see text S2). 

In contrast, the relatively blue-shifted structures of the Na I D lines show a rather complex 
evolution. The number of features, their intensity and width are difficult to estabUsh. Nev- 
ertheless, for the sake of discussion, four main components, which we will indicate as A, B, 
C and D, can be tentatively identified in the first two epochs (Fig. 2). Components B, C and 
D strengthen between day —2 and day +14 while component A remains constant during this 
time interval. The situation becomes more complicated on day +61: components C and D 
clearly start to decrease in strength, but component B remains almost constant and component 
A becomes definitely deeper and is accompanied by a wide absorption that extends down to a 
rest- frame heliocentric velocity ~ —50 km s~^ (Fig. 1, Fig. S2). After this epoch there is no 
evidence of evolution, and component A remains the most intense feature up to the last phase 
covered by our observations, more than four months after the explosion. 

Variable interstellar absorption on comparably short timescales has been claimed for some 
Gamma Ray Bursts (ORB), and it has been attributed by some authors to line-of-sight geomet- 
rical effects, due to the fast GRB expansion coupled to the patchy nature of the intervening 
absorbing clouds (18). Our data clearly show that despite the marked evolution in the Na I D 
lines, Ca II H&K components do not change with time (see Fig. 2, Fig. S3 and the discussion in 
sections S3 and S4). Therefore, in the case of SN 2006X, transverse motions in the absorbing 
material and line-of-sight effects due to the fast SN photosphere expansion (typically 10^ km 
s~^) can be definitely excluded, since they would cause variations in all absorption features. 

For this reason we conclude that the Na I features seen in SN 2006X, arising in a number 
of expanding shells (or clumps), evolve because of changes in the CSM ionization conditions 
induced by the variable SN radiation field. In this context, the different behavior seen in the 
Na I and Ca II lines is explained in terms of the lower ionization potential of Na I (5.1 eV, 
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corresponding to 2417A) with respect to Ca II (11.9 eV, corresponding to 1045 A), their dif- 
ferent recombination coefficients and photoionization cross sections coupled to a UV-deficient 
radiation field (see text S4). Regrettably, not much is known about the UV emission of SNe 
la shortwards of 11 00 A (8, 19). From a theoretical point of view, a severe UV line blocking 
by heavy elements like Fe, Co and Mg is expected (20). An estimate of the Na 1 ionizing flux, 
SuV', can be derived from a synthetic spectrum of a Type la SN at maximum light (21 ), and this 
turns out to be Suv ~5x 10^° photons s~^. One can verify that this flux is largely sufficient to 
fully ionize Na 1 up to rather large distances ~5x 10^^ cm. 

Nevertheless, since the recombination timescale must be of the order of 10 days, this 
requires an electron density Ue as large as 10^ cm~^ ( see section S4). Given the low abundance 
of any other element, such a high electron density can be produced only by partial hydrogen 
ionization. Due to the severe line blocking suffered by Type la SNe (20), the flux of photons 
capable of ionizing H is very small (~4x 10^^ photons s~^) and this imposes that the gas where 
the Na I time-dependent absorptions arise must be confined within a few 10^^ cm from the SN 
(see text S4). In a SN of this t5^e, the flux in the 11 20-2640 A band decreases by a factor 
of ten in the first two weeks after maximum light (8, 19). Since at a distance of ~10^^ cm 
from the SN the ionization timescale Tj for Na I is much shorter than Tr, the ionization fraction 
grows with time following the increase of the UV flux during the pre-maximum phase, while 
after maximum it decreases following the recombination timescale. This would explain the 
overall growth of the blue components' depth, as shown by our data, in terms of an increasing 
fraction of neutral Na, while the different evolution of individual components would be dictated 
by differences in the densities and distances from the SN. Moreover, once all the Na II has 
recombined (which should happen within a few r^, i.e. ~1 month), there should be no further 
evolution, in qualitative agreement with the observations. Additionally, since the flux of photons 
that can ionize Ca n is more than four orders of magnitude less than in the case of Na I (see 
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section S4), the corresponding ionization fraction is expected to be of a few per cent only. 
Therefore, the recombination of Ca III to Ca II does not produce measurable effects on the 
depth of Ca II H&K lines, as it is indeed observed (see text S3). 

The H mass (M(H)) contained in the shells generating the observed absorptions can be es- 
timated from our observations after making some conservative assumptions. The Na I column 
density A^(Nal) deduced from the most intense feature (D, day +14) is A^(Nal) ~ 10^^ cm~^. 
Assuming that the material generating this component is homogeneously distributed in a thin 
spherical shell with radius 10^^ cm, a solar Na/H ratio (log Na/H=— 6.3) and complete Na re- 
combination, an upper limit to the shell mass can be estimated as M(H) < 3 x 10~^ Mq (this 
value is reduced by a factor 100 for material at about 10^^ cm, the most likely distance for com- 
ponents C and D; see below). Even in the case of complete ionization, such a H mass would 
produce an Ha luminosity of ~4xl0^'' erg s~^, which is two orders of magnitude below the 
3-sigma upper limits set by our observations at all epochs (Table S2) and by any other SN la 
observed so far (6). Therefore, the absence of narrow emission lines above the detection limit 
is not in contradiction with the presence of partially ionized H up to masses of the order of 0.01 
Mo. 

However, photo-ionization alone cannot account for the fact that not all features increase 
in depth with time (Fig. 2). In fact, on day +61, components C and D turn back to the same 
low intensity they had on day —2. One possible explanation is that the gas is re-ionized by 
some other mechanism, like the ejecta-CSM interaction. In this case, the absorbing material 
generating components C and D must be close enough to the SN so that the ejecta can reach it 
in about one month after the explosion (~10^^ cm for maximum ejecta velocities of 4x 10^ km 
s~^). Similarly, in order not to be reached by the ejecta more than four months after. A, B and 
the broad high-velocity components must arise at larger distances (>5 x 10^^ cm). This scenario 
is not ruled out by the lack of radio emission from SN 2006X (15). In fact, in the light of our 
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Figure 2: Evolution of the Na I D2 and Ca II K line profiles between day —2 (black), day +14 
(red) and day +61 (blue, Na I D2 only). The vertical dotted lines mark the four main variable 
components at —3 (A), +20 (B), +38 (C) and +45 (D) km s^^. For comparison, the upper panel 
shows the R{0), R{1) and P(l) line profiles of the (0-0) vibrational band of the CN E^S -X^S. 
The velocity scale refers to the R{0) transition (3874.608 A). Color coding is as for the other 
two panels. 
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current understanding of the ejecta-CSM interaction mechanism (22), the presence of similar 
shells with masses smaller than a few 10^'^ M©, cannot be excluded by radio non-detections of 
SNe la in general (7). Our findings are consistent with upper limits on the radio flux set by our 
VLA observations, obtained about 10 months after the explosion (see section SI), which are 
comparable to the best upper limits set on the radio luminosity of other normal SNe la (7). 

If we adopt the velocity of the CN lines as indicative of the host galaxy rotation component 
along the line of sight at the SN location, then our observations provide soUd evidence of CSM 
expanding at velocities that span a range of about 100 km s~^ (Fig. 2). 

The most important implication of these observations is that they show that this circumstellar 
material was ejected from the progenitor system in the recent past. For instance, with a shell 
radius of 10^^ cm and a wind velocity of ~50 km s~^ , the material would have been ejected some 
50 years before the explosion. This almost certainly rules out a double-degenerate scenario for 
SN 2006X, where the supernova would have been triggered by the merger of two CO white 
dwarfs. In this case, no significant mass loss would be expected in the phase immediately 
preceding the supernova. Thus, a single-degenerate model is the favored model for SN 2006X, 
where the progenitor accreted from a non-degenerate companion star. 

Mean velocities for the circumstellar material of ~50 km s~^ are comparable to those re- 
ported for the winds of early red giant (RG) stars (22); velocities matching our observations 
are also expected for late subgiants. The observed material is moving more slowly than would 
be expected for winds from main sequence donor stars or from compact helium stars. Of the 
two major formation channels proposed for SN la with a non-degenerate donor star (23), these 
wind velocities seem more consistent with the shorter-period end of the "symbiotic" formation 
channel. The observed structure of the circumstellar material could be due to variability in the 
wind from the companion RG; considerable variability of RG mass loss is generally expected 
(24). 
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An alternative interpretation of these distinct features is that they arise in the remnant shells 
of successive novae, which can create dense shells in the slow moving material released by the 
companion (25, 26). This seems to require an aspherical shell geometry in order to match the 
observed low velocities (see S6 for further details). Not only might this be expected a priori 
(27), observations of the 2006 outburst of RS Ophiuchi show that there is an equatorial density 
enhancement which strongly restrains the expansion of the nova shell (28, 29, 30). 

One crucial issue is whether what we have seen in SN 2006X represents the rule or is rather 
an exceptional case. Other cases of SNe la showing negative velocity components are known, 
like SN 199 IT and SN 1998es (see Fig. S5 and the discussion in S5). Unfortunately, multi- 
epoch high-resolution spectroscopy is not available for these objects (to our knowledge, the 
SN 2006X data set is unique in this respect), and therefore time variability cannot be demon- 
strated. Nevertheless, the data clearly show components approaching the observer at velocities 
which reach at least 50 km s~^ with respect to the deep absorption that we infer to be produced 
within the disks of the respective host galaxies. This, and the fact, that SN 2006X has shown no 
optical, UV and radio peculiarity whatsoever, supports the conclusion that what we have wit- 
nessed for this object is common to normal SN la, possibly all of them, even though variations 
due to different inclinations of the line of sight with respect to the orbital plane may exist. 
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SI. Materials and Methods 

We have observed SN 2006X on 5 epochs spanning about 4 months (Table SI). High res- 
olution spectra were obtained with the European Southern Observatory's (ESO) Very large 
Telescope (VLT) on Cerro Paranal (Chile), equipped with the Ultraviolet and Visual Echelle 
Spectrograph (UVES) (SI), and with the Keck I Telescope equipped with the High Resolution 
Echelle Spectrometer (HIRES) (S2). For UVES we have used the 390+580 setting, which cov- 
ers simultaneously three wavelength ranges (3290-4500A, 4780-5740A and 5830-6800A), with 
a full width half maximum (FWHM) resolution of 7 km s~^. In the case of HIRES, we have 
used a setting optimized for the Nal D region, that covers the wavelength range 3900-8350A 
with a FWHM resolution of 6.8 km s^. UVES data have been reduced using the UVES Data 
Reduction Pipeline (S3), while HIRES data have been processed using standard procedures for 
Echelle spectra. Wavelength calibration has been achieved using Thorium- Argon lamps. The 
final RMS accuracy is about 0.15 km s~^. The wavelength scale was corrected to the rest-frame 
adopting a host galaxy recession velocity of 1571 km s^^ (S4). To compensate for the Earth's 
motion, a heliocentric velocity correction has been applied to the data (Table SI). The effect of 
atmospheric lines has been checked using a spectroscopically featureless bright star (HR 3239) 
observed with the same instrumental setup as for the science data. The spectral region of inter- 
est turns out to be essentially free of telluric features. Therefore, the broad absorption visible 
starting with day +61 and reaching Vh ~ —70 km s~^ is real and not affected by atmospheric 
lines. Na I column densities have been estimated fitting Voigt line profiles with VPFIT (S5). 
Finally, epochs have been computed from B maximum light, which took place on February 20, 
2006 (S6). 

In order to set upper limits to Ha and He 1 5876A luminosities (Table S2), the UVES spectra 
have been calibrated by means of a reference response function. Since the observations have 
been obtained under sky transparency conditions that ranged from clear/photometric (epochs 
—2, +14, and +121) to thin cirrus (+61), taking into account the high instrument stabiUty, the 
expected flux calibration accuracy is of the order of 20-30%. To the best of our knowledge, the 
upper limits on Ha and He I 5876A presented here are the latest ever published. The signal-to- 
noise ratio on the continuum in the Na I D Unes region ranges from ~70 (day —2) to ~20 (day 
+121). 

To confirm the low mass of CSM estimated from the Na I observations, we have observed 
SN 2006X with the Very Large Array ( S8), in the C configuration, at 6 cm (4.8 GHz) on Novem- 
ber 17 and at 3.6 cm (8.4 GHz) on November 20, 2006, corresponding to +270 and +273 days 
after B maximum light, respectively. The SN was not detected at both wavelengths, and the 2 
SD limits are 0.07 mJy and 0.09 mJy for the two bands respectively (S9), which are comparable 
to the best limits available so far (570). At the epoch sampled by our observations, the SN ejecta 
have reached a distance of ~10^^ cm from the explosion site. Our upper limits on the radio flux 
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UT Date 
(2006) 


Phase 
(days) 


Instrument/Telescope 


Total Integration Time 
(seconds) 


Heliocentric Correction 

(km s"^) 


18/02 


-2 


UVES/VLT 


4175 


+14.6 


06/03 


+14 


UVES/VLT 


8940 


+7.2 


22/04 


+61 


UVES/VLT 


15025 


-15.4 


06/06 


+105 


HIRES/KECK 


3600 


+5.0 


25/06 


+121 


UVES/VLT 


15025 


-28.0 



Table SI: High resolution spectroscopic observations of SN 2006X. Phase refers to S-band 
maximum Ught, attained on February 20, 2006 (S6) 



constrain the corresponding CSM mass to be smaller than a few 10^^ solar masses (M0), which 
is fully consistent with the estimates deduced from the observed Na 1 absorptions. Very similar 
upper limits have been set for standard SN la events like SN 198 IB, 1989M, 1998bu and 1992A 
(SIO). 

From the low resolution spectroscopic data that we have obtained in parallel to the high 
resolution data set discussed in this paper, we confirm the results of the early classification 
spectroscopy (SI 1). SN 2006X appears to be a normal la, very similar to SN 2002bo (S12), with 
somewhat higher photospheric velocities. The derivation of photometric parameters is hindered 
by a rather strong reddening (E{B — V) >1.1) and an anomalous extinction law (Ry <2). 
However, the decline rate, the light and color curve shapes and the absolute magnitude are 
within the range of normal Type la SNe (S6). 

S2. The interstellar material 

The apparent position of SN 2006X is very close to a spiral arm of NGC 4321, which is 
inclined by 28 degrees with respect to the line of sight (SI 3). Due to the relatively high extinc- 
tion (Sll), the SN most likely exploded within or behind the disk of the host galaxy. Moreover, 
the SN is projected onto the receding side of the galaxy, and the component of the rotation 
velocity along the line of sight at the apparent SN location is about 75 km s~^ (S4, S13), which 
coincides with the strongly saturated Na 1 D component (Fig. SI). This suggests rather un- 
equivocally that the time-invariant, very deep absorptions arise within the disk of NGC 4321. 
The equivalent widths of the Na 1 D lines, measured in the first epoch spectrum, are 670±5 and 
625 ±4 mA for the D2 and Di components respectively. The implied total Na I column density, 
estimated fitting multiple Voigt line profiles with VPFIT (S5), is as large as log A?^ ~14.3. For 
a Milky Way-like dust mixture this would turn into a color excess Eb-v ~1-1 (S14). More- 
over, the CN lines clearly visible in our first two spectra at t'/i=73.6±1.0 km s^^ (weighted 
average of the two epochs) have an unprecedented depth (see also ref. S15). On the first epoch 
spectrum the equivalent widths for the three lines are: EW[R{0)]=91±2, EW[R{1)]=60±3 
and £'1^[P(1)]=42±3 mA that, after correcting for saturation effects, correspond to column 
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Phase Ha He 1 5876 





7 km s ^ 


50kms-^ 


7 km s ^ 


50kms-^ 


-2 


2.2 


16.0 


2.4 


17.0 


+14 


1.1 


8.0 


1.1 


7.6 


+61 


0.6 


4.0 


0.6 


4.4 


+121 


0.3 


2.2 


0.4 


2.6 



Table S2: Upper limits (3 SD) for Ha and He I 5876A luminosities measured on the UVES 
spectra. Luminosities are expressed in 10^^ erg s"^ and were computed for two Gaussian line 
profiles with FWHM 7 km ~^ (matching the instrumental resolution) and 50 km s"^^, assuming 
a distance of 16.1 Mpc for NGC 4321 (S7). 

densities of log[iV^(o)]=14.0±0.1, log[iV/j(i)]=13.6±0.1 and log[A^R(i)]=13.6±0.1 respectively. 
Stars having comparable column densities in our own Galaxy are known to have color excesses 
larger than Eb-v=^ (S16). Applying the usual curve of growth method (S16), an excitation 
temperature of 3.0 ± 0.2 °K is derived for this cloud. 

A closer inspection of Fig. SI shows that the deep Na 1 and Ca II absorptions are not due 
to a single cloud. In fact the Ca 11 H line, which is less affected by saturation, clearly displays 
at least two components that differ in velocity by ~10 km s"^. On the contrary, the CN lines 
appear as single, unresolved features (the velocity dispersion parameter is b ~2.3 km s~^), 
indicating that this molecular absorption arises within an inner subset (as seen along the line of 
sight) of multiple clouds that most likely contribute to the strong Na 1 and Ca 11 absorptions. 

The peculiarity of the interstellar material is indicated also by the unusual wavelength de- 
pendency of the continuum linear polarization shown by SN 2006X. This is in fact significantly 
different from that typically observed in extinguished stars within our own Galaxy (SI 8). 

The presence of strongly saturated Na ID interstellar components could in principle mask 
the existence of weaker features, which happen to fall in the same velocity range but arise in 
a completely different region. Therefore, the presence of an absorption feature at the typical 
velocity of a RG wind (10 km s~^) cannot be ruled out. Only the study of non-extinguished 
objects, i.e. free of strong interstellar absorptions, will allow us to detect slow moving material 
in the surroundings of Type la SNe. 

S3. Lack of evolution in Ca II H&K lines 

The Nal D evolution is very clearly displayed by both Di and D2 lines, which show prac- 
tically identical profiles at all epochs (Fig. S2). SUght saturation effects are visible in the most 
intense components of the D2 line. 

The behavior of the Ca II H & K lines is radically different from that of the Na I D lines. 
On the first epoch (day —2), the profiles of Na I D2 and Ca n K show very similar intensity 
and velocity profiles (Fig. S3). In particular, components A and B can be easily identified. The 
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Figure SI: Comparison between the profiles of CN R(0), Na I Di and Ca II K lines for the 
first UVES epoch (day —2). The vertical dashed line marks the velocity of the R(0) absorption 
center {vh='13.6 km s"^). 
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differences become dramatic at the second epoch (day +14), when Ca II K remains practically 
unchanged while Na 1 D2 develops deep absorptions for components B, C and D. The signal- 
to-noise ratio in the Ca 11 region is lower, but a variation similar to that taking place for the Na 1 
lines would have been definitely detected, at least for components B, C and D. Due to the SN 
fading and its spectral evolution, the continuum in the Ca II H&K region becomes very weak as 
time goes by, making the detection of the K line more and more difficult. Even heavily binning 
our third epoch data (day +61), the signal-to-noise ratio achieved on the resulting profile is 
not sufficient to draw any firm conclusion (Fig. S3). At later epochs, due to the very low SN 
continuum, the H&K lines could not be detected. 

Geometric effects related to the fast expansion of the radiation ring and the consequent 
change in the covering factor coupled to the patchy structure of the intervening absorbing 
clouds, have been advocated to explain the time variability claimed for some Gamma Ray Bursts 
(GRB) (S19, S20). What emerges from this analysis is that a Ca 11 H&K evolution at the level 
seen in the Na I D lines can be definitely excluded. In turn, this rules out transverse motions or 
line-of-sight effects as responsibles for the evolution observed in SN 2006X. 



S4. The location of the absorbing material in SN 2006X 

After sodium is ionized by the SN radiation field, which drops quickly after maximum light 
(S21), it starts recombining with a time scale which depends on the electron number density 
He and the radiative recombination coefficient /3: Tr = {ue Given the time scale of 

the observed changes in the NalD features, we estimate r,. ~10 days, which implies ~ 10^ 
cm~'^ for (3 ~5 x 1 0^ cm"^ s~ ^ ( S22 ). Such an high electron density can only come from partial 
ionization of hydrogen. Therefore, the maximum distance for the gas where the time-variable 
Nal features arise is dictated by the ability of the SN to ionize hydrogen, i.e. by the flux Suv of 
photons at wavelengths shorter than 912A, times the duration AtgN of the UV emitting phase. 

An upper limit to the number of hydrogen atoms that are ionized is therefore given by the 
total number of ionizing photons, Suv^tsN, which in turn must be equal to Vh Ue, with Vh 
being the volume of the layer responsible for the varying absorption. Thus, an upper limit to 
this volume is given by: 

Vh < ^MX^, (51) 

Tie 

Assuming that the absorbing material is confined within a thin shell of thickness Ar (with 
Ar <^ Vh), the maximum radius of such a shell is: 



Suv ^tsN 



1/3 

(32) 



Att He Ar/r 

Using Suv= 4.4x 10'^^ photons s^^, the value derived from the synthetic spectrum at max- 
imum light (S23), and a light curve width of AtsN= 20 days, one gets rn <4xl0^^ cm for 
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Figure S2: Comparison between the Nal Di and D2 line profiles on days —2, +14, +61 and 
+121. For presentation, the intensity of the Di line has been multiplied by 2.0, that is the D2/D1 
ratio expected from the spin orbit statistical weights of the Na(^P3/2) and Na(^Pi/2) transitions. 
In each panel, the dotted curve traces the atmospheric absorption spectrum in the D2 region. 
In order to account for the different heliocentric corrections applied to the SN data, the telluric 
absorption spectrum has been shifted in velocity to match the atmospheric features visible in 
the SN spectra. The vertical dashed lines mark the main velocity components. 
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Restframe Heliocentric Velocity v,^ (km s"^) 

Figure S3: Comparison between the evolution of NaT D2 and Ca II K lines during the first three 
epochs covered by our observations. In order to increase the signal-to-noise ratio on day -1-6 1 
the Ca II K data have been binned to a 0.3 A step, corresponding to about 23 km s^. The error 
bars indicate ±1 SD level, estimated from the data within each single bin, composed by 21 
original data points. In each panel, the dotted curve traces the atmospheric absorption spectrum 
in the D2 region. 
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ne=lO^ cm~^ and Ar/r=0.01. The largest uncertainty in estimating rn comes from the total 
fluence of ionizing photons (Suv^tsN)- In fact, the synthetic spectrum we have used is a best 
fit to an HST spectrum of SN 1992A extending down to about I6OOA and obtained 6 days past 
maximum light and (S24), i.e. when the UV flux has significantly dropped (S21 ). Nevertheless, 
thanks to the cubic root dependency on it, even a 10^ times higher fluence would result in an 
upper limit of ~5xlO^^ cm. This clearly demonstrates that the material responsible for the 
time-variable features is confined well within the circum- stellar domain. 

This creates a very marked distinction between Type la SNe and GRBs, for some of which 
time-dependent UV absorption features have been claimed (S19, S20, S25). In fact, GRBs 
and their early after-glows have a very strong X-ray/UV radiation field, so that all possible 
circumstellar gas is completely ionized. For instance, GRB 021004 was shown to be able to 
completely ionize the interstellar material out to about 100 pc (~3x 10^° cm) (S19), while for 
GRB 060418 the ionization was found to reach ~1.7 kpc (~5x 10^^ cm) (S25). 

As far as calcium is concerned, we note that the flux of photons capable of ionizing Ca II 
computed from the synthetic spectrum (S23) is ~3.5x 10^^ photons s~^. This is more than four 
orders of magnitude smaller than the corresponding Na I ionizing flux. As a consequence, the 
ionization fraction of Ca II is significantly smaller and hence, as it recombines, no detectable 
effect on the depth of Ca II H&K lines is expected, as shown by the observations (see S3 and 
Fig. S3). The same mechanism has been proposed to explain the inter- stellar line variability 
observed in galactic stars which, due to the much lower densities, takes place on timescales of 
years to decades (S26). 

An alternative explanation for the presence of expanding shells in the immediate surround- 
ings of SN 2006X is that they have been lost by a nearby star, physically not related to the 
progenitor system but close enough to be influenced by the SN radiation field. In fact, shells 
expanding with velocities up to several tens of km s~^ are known to be present, for instance, 
around many galactic OB associations (S26). Since the material has to be relatively close to the 
SN in order to be influenced by its UV radiation, this alternative scenario requires that the SN 
exploded within a cluster or a star association. 

In order to study the immediate surroundings of the explosion site, we have used a deep pre- 
explosion image obtained with the ESO-VLT FOcal Reducer and low dispersion Spectrograph 
(FORSl), (S28). Observations were performed on May 12, 1999 in the R passband (S28), 
with an exposure time of 600 seconds and a seeing of 0.8 arcsec (FWHM), corresponding to 
a limiting magnitude R ~24.0 (5 SD level) at the SN location (Fig. S4, left panel). To better 
characterize the SN environment, we have searched the Hubble Space Telescope Archive for 
high resolution imaging. The only suitable data we could find were obtained with the Advanced 
Camera for Surveys (ACS), (S29) on May 21, 2006 (S30) in high resolution mode (0.027 arcsec 
pixer^) with the F435W (1480 seconds), F555W and F775W (1080 seconds) filters. In order 
to increase the signal-to-noise ratio, the three images have been registered and stacked (Fig. S4, 
right panel). The resulting image reveals a wealth of resolved objects, which can be identified 
as red, blue super-giants, groups of O stars and OB associations, with typical sizes of ~100 pc. 
All such associations and stellar groups (like those marked with labels A and B on the ACS 
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Figure S4: The host galaxy imaged before (left, VLT-FORSl) and after (right, HST-ACS) the 
explosion of SN 2006X. The field of view is 26.5 x 22.5 arcsec for both images. The SN position 
in the FORSl image is marked by a circle. The spatial scale on the ACS images is about 2.1 pc 
pixel" ^. 

image) are clearly detected in the deep FORSl image, which shows also some of the brightest 
single stars. 

Unfortunately, the ACS image was obtained after the explosion, when the SN was still 
bright, and hence it is not possible to directly inspect the SN surroundings. Nevertheless, the 
presence of groups like that marked with A in Fig. S4 can be excluded. This is confirmed by 
the deep FORSl pre-explosion image, which does not show any trace of stellar associations 
(like A) or even more compact groups (like B) at the SN location. This makes the presence 
of an isolated star losing material in the immediate vicinity of the SN progenitor system very 
unlikely. 

S5. Evidences of blue-shifted Na I absorption lines in other SNe la 

Due to the large required exposure times, high resolution spectroscopy of SN la is rather 
rare. Additionally, since the main aim of those investigations was the study of the supposedly 
time-invariant inter-stellar medium, data were usually obtained on one single epoch around 
maximum light (S31, S32). Only in recent times, with the advent of 8- 10m class telescopes, 
multi-epoch high resolution spectroscopy has been performed, with the aim of detecting possi- 
ble emission lines produced by the SN ejecta-CSM interaction (S33, S34). Nevertheless, due to 
the underlying scientific driver of these projects, observations were carried out during the pre- 
and near-maximum phase, when the Na I D features are at their minimum intensity in SN 2006X 
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(Fig. SI). For example, the spectra of SN 2001el, observed some days before maximum light, 
do not show any clear signs of Na 1 D blue components. Interestingly, they do show at least two 
blue components in the Ca II H profile at ~ — 18 and —34 km s~^ with respect to the deep Na I 
D absorption (S3 3). 

In the course of our long term nearby SN monitoring, we had the chance to obtain high 
resolution spectroscopy of two SN la, namely SN 1991T (S35) and SN 1998es (S36). Both 
spectra are unpublished. SN 199 IT was observed on Jun 08, 1991 (day +40) with the ESO 
1.4m Coude Auxiliary Telescope (CAT) equipped with the Coude Echelle Spectrograph (S37), 
while the spectrum of SN 1998es was obtained on Nov 24, 1998 (day -1) with the ESO 1.5m 
telescope equipped with the Fiber-fed Extended Range Optical Spectrograph ( S3 8). The Na I D 
line profiles of these two objects are remarkably similar to that of SN 2006X (Fig. S5). Besides 
a strong absorption, which is most likely generated within the disk of the host galaxy, some blue 
components are clearly visible. 

To the best of our knowledge, the only other SN for which a data set somewhat similar 
to that of SN 2006X has been obtained is the core-collapse. Type Iln SN 1998S (S39). High 
resolution spectroscopy of this object has revealed a number of Na I D components within the 
host galaxy, the bluest of which (vh ~ —100 km s~^) deepened significantly during the 19 
days spanned by the two epochs available, implying a Na 1 column density increase of about 
1 dex. This fact, together with the detection of time-evolving narrow H and He P-Cyg profiles 
at the same Na I velocity, was interpreted as a signature of the outflows from the super-giant 
progenitor of SN 1998S, arising in a dense shell, expanding at about 50 km s~^ (S39, S40). 



S6. The CSM structure around SN 2006X: swept-up nova shells? 

A possible interpretation of the distinct features in the CSM is that they arise in the remnant 

shells of successive novae (S41). The recurrent ejection of a few 10~^ Mq of H at velocities of 
~4000 km s~^ can create dense shells in the slow moving material released by the companion, 
evacuating significant volumes around the progenitor system (S42). This explanation requires 
that the high-velocity nova ejecta have, by sweeping up the stellar wind of the donor star, been 
slowed down to speeds several times lower than the asymptotic shell velocities seen, for in- 
stance, in RS Ophiuchi, i.e. ~ 300 km s"^ (S43, S44). For SN 2006X, if we assume that the 
observed shells have all been swept up in the energy-conserving phase (S45), and if we take an 
upper limit of the wind-mass loss rate from the companion, Mwind> as 10"^ Mq yr~^ from the 
radio observations, we obtain an upper limit on the mass ejected in each nova outburst: 

M „ T r>— 8 71 ^ -^^wind f ^hell \ f ^^ova \ ^^nova /oo\ 
nova < 2 X 10 Mq—— — 31 n -{7, — . (S3) 

10 'M0yr ^ V75kms V V 1000 km s / 10 yr 

where M and V refer to mass and velocity, and the subscripts nova and shell refer to the mass 
ejected in each nova outburst and the shell formed as the nova ejecta sweep up the wind, whilst 
Atnova represents the nova recurrence time. 
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Figure S5: Comparison between the Sodium D2 line profiles for three SN la: 1991T (upper), 
1998es (middle) and 2006X (lower). For the sake of clarity, spectra have been arbitrarily shifted 
in the horizontal direction in order to place the center of the strongest D2 component at zero 
velocity. 
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Although this value is around an order of magnitude smaller than might be expected from re- 
cent nova calculations, we note that the ejecta mass is a steeply decreasing quantity as the white 
dwarf approaches explosion (S41, S46). In addition, the estimate above assumes that the wind 
from the companion is spherically symmetric. If the wind is concentrated towards the orbital 
plane, and if we are observing the shell structure close to the orbital plane, the shells in our Une- 
of-sight would be slowed down more than in the spherically symmetric case. The discrepancy 
would be further reduced if the nova shell is mainly slowed down in the momentum-conserving 
phase rather than the energy-conserving phase, which is made more likely in this asymmet- 
ric case where the mass is concentrated in the orbital plane. Exactly such an axisymmetric 
geometry and rapid deceleration is suggested by recent observations of RS Oph (S47). 

Indeed, the fact that there is neutral Na around the progenitor of SN 2006X requires cool ma- 
terial and suggests that the shell must have cooled down sufficiently and be in the momentum- 
conserving phase. 

If the recurrent nova interpretation requires an aspherical geometry, then this implies that 
the apparent shell structure with the velocities we have observed would only be seen in SNe 
la which are observed close to the orbital plane of the system. This may be a useful future 
discriminant for this model. We also suggest that this interpretation is made more satisfying 
if the loss of the strong Na line at ~40 km s^^ (components C and D in Fig. 2) is due to the 
interaction of the supernova ejecta with that nova shell, as that places a useful scale on the 
distance of the shell from the supernova. As a typical SN velocity of ~10^ km s~^ is some 250 
times faster than the shell, the interaction after 1 month corresponds to a recurrence time of ~25 
yr, which is a plausible value. The actual recurrence time could well be less, as that shell is not 
guaranteed to be the innermost shell, merely the innermost visible one. 

As a check on the consistency of this model, we can combine our inferred electron density 
of He ~ 10^ cm~^ with the dimensions of the nova shell to produce an estimate of the mass of 
the shell: 

(^) ~ - X (t) (lo^) (xo^) i^f' 

where, somewhat arbitrarily, we have taken a shell thickness (Ar) that is ~ 1% of the radius 
of the shell. Our greatest uncertainty is probably contained within the ionisation fraction X, 
which is needed to obtain the mass per free electron, though the quadratic dependence on rsheii 
alone leads to considerable uncertainty in the estimate of Mghcu- Note that we have assumed the 
shell is spherical; a wind concentrated in the equatorial plane may be ~ 10 times less massive. 

If we assume that the mass in the shell is dominated by mass swept up from the wind (i.e. 
Mwind ^ Mnova) ^nd then take the shell radius to be the product of the recurrence time and shell 
velocity we can write: 

[Mq/jt) ~ ^ \xj VlOi3cmy Vl05cm-V 1^50 km/s/ 1^ 10 yr j ' 

The wind loss rate derived in this way is consistent with our expectation in Equation S3. 
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